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Abstract — This paper presents a novel adaptive reduced-rank 
multi-input multi-output (MIMO) equalization scheme and algo- 
rithms based on alternating optimization design techniques for 
MIMO spatial multiplexing systems. The proposed reduced-rank 
equalization structure consists of a joint iterative optimization 
of two equalization stages, namely, a transformation matrix that 
performs dimensionality reduction and a reduced-rank estimator 
that retrieves the desired transmitted symbol. The proposed 
reduced-rank architecture is incorporated into an equalization 
structure that allows both decision feedback and linear schemes 
for mitigating the inter-antenna and inter-symbol interference. 
We develop alternating least squares (LS) expressions for the de- 
sign of the transformation matrix and the reduced-rank estimator 
along with computationally efficient alternating recursive least 
squares (RLS) adaptive estimation algorithms. We then present 
an algorithm for automatically adjusting the model order of the 
proposed scheme. An analysis of the LS algorithms is carried 
out along with sufficient conditions for convergence and a proof 
of convergence of the proposed algorithms to the reduced-rank 
Wiener filter. Simulations show that the proposed equalization 
algorithms outperform the existing reduced-rank and full-rank 
algorithms, while requiring a comparable computational cost. 

Index Terms — MIMO systems, equalization structures, param- 
eter estimation, reduced-rank schemes. 



I. Introduction 

THE high demand for performance and capacity in wire- 
less networks has led to the development of numerous 
signal processing and communications techniques for employ- 
ing the resources efficiently. Recent results on information 
theory have shown that it is possible to achieve high spectral 
efficiency fT| and to make wireless links more reliable J3], 
J5) through the deployment of multiple antennas at both 
transmitter and receiver. In multi-input multi-output (MIMO) 
communications systems, the received signal is composed 
by the sum of several transmitted signals which share the 
propagation environment and are subject to multi-path propa- 
gation effects and noise at the receiver. The multipath channel 
originates inter-symbol interference (ISI), whereas the non- 
orthogonality among the signals transmitted gives rise to inter- 
antennas interference (IAI) at the receiver. 
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In order to mitigate the effects of ISI and IAI that reduce 
the performance and the capacity of MIMO systems the 
designer has to construct a MIMO equalizer. The optimal 
MIMO equalizer known as the maximum likelihood sequence 
estimation (MLSE) receiver was originally developed in the 
context of multiuser detection in (4). However, the exponential 
complexity of the optimal MIMO equalizer makes its imple- 
mentation costly for multipath channels with severe ISI and 
MIMO systems with many antennas. In practice, designers 
often prefer the deployment of low-complexity MIMO re- 
ceivers such as the linear |5|. |6'. the successive interference 
cancellation-based VBLAST and decision feedback equal- 
izers (DFE) Il8l-lfl4l The DFE schemes HD-CH) can achieve 
significantly better performance than linear ones due to the 
interference cancellation capabilities of the feedback section. 
These receivers require the estimation of the coefficients used 
for combining the received data and extracting the desired 
transmitted symbols. A challenging problem in MIMO systems 
lfT31l is encountered when the number of elements in the 
equalizer or the number of antenna pairs is large, which 
is key to future applications |[T6l - lfT8l . In these situations, 
an estimation algorithm requires substantial training for the 
MIMO equalizer and a large number of received symbols to 
reach its steady-state behavior. 

There are many algorithms for designing MIMO equalizers, 
which possess different trade-offs between performance and 
complexity lfl9l . In this regard, least squares (LS)-based 
algorithms are often the preferred choice with respect to con- 
vergence performance. However, when the number of filter el- 
ements in the equalizer is large, an adaptive LS-type algorithm 
requires a large number of samples to reach its steady-state 
behavior and may encounter problems in tracking the desired 
signal. Reduced-rank techniques 11201 - 1341 are powerful and 
effective approaches in low-sample support situations and in 
problems with large filters. These algorithms can exploit the 
low-rank nature of signals that are found in MIMO commu- 
nications 136] in order to achieve faster convergence speed, 
increased robustness against interference and better tracking 
performance than full-rank techniques. By projecting the input 
data onto a low-rank subspace associated with the signals of 
interest, reduced-rank methods can eliminate the interference 
that lies in the noise subspace and perform denoising ||201 - ||34| . 
Prior work on reduced-rank estimators for MIMO systems is 
extremely limited and relatively unexplored, being the work 
of Sun et al. |25ll one of the few existing ones in the area. A 
comprehensive study of reduced-rank equalization algorithms 
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for MIMO systems has not been considered so far. It is 
well known that the optimal reduced-rank approach is based 
on the eigen-decomposition (EVD) of the known input data 
covariance matrix R l20l . However, this co variance matrix 
must be estimated. The approach taken to estimate R and 
perform dimensionality reduction is of central importance and 
plays a key role in the performance of the system. Numerous 
reduced-rank strategies have been proposed in the last two 
decades. The first methods were based on the EVD of time- 
averaged estimates of R 1201 . in which the dimensionality 
reduction is carried out by a transformation matrix formed by 
appropriately selected eigenvectors computed with the EVD. 
A more recent and elegant approach to the problem was taken 
with the advent of the multistage Wiener filter (MSWF) 11221 . 
which was later extended to adaptive versions in 11231 . Il24l . 
and MIMO applications (25l . Another related method is the 
auxiliary vector filtering (AVF) algorithm Il26l - ll28l . which can 
outperform the MSWF. A key limitation with prior art is the 
deficient exchange of information between the dimensionality 
reduction task and the subsequent reduced-rank estimation. 

In this work, we propose adaptive reduced-rank MIMO 
equalization algorithms based on alternating optimization de- 
sign techniques for MIMO spatial multiplexing systems. The 
proposed reduced-rank equalization structure and algorithms 
consist of a joint iterative optimization that alternates between 
two equalization stages, namely, a transformation matrix that 
performs dimensionality reduction and a reduced-rank esti- 
mator that suppresses the IAI caused by the associated data 
streams and retrieves the desired transmitted symbol. The 
essence of the proposed scheme is to change the role of the 
equalization filters and promote the exchange of information 
between the dimensionality reduction and the reduced-rank 
estimation tasks in an alternated way. In order to estimate the 
coefficients of the proposed MIMO reduced-rank equalizers, 
we develop alternating least squares (LS) optimization algo- 
rithms and expressions for the joint design of the transforma- 
tion matrix and the reduced-rank filter. We derive alternating 
recursive LS (RLS) adaptive algorithms for their computation- 
ally efficient implementation and present a complexity study 
of the proposed and existing algorithms. We also describe 
an algorithm for automatically adjusting the model order of 
the proposed reduced-rank MIMO equalization schemes. An 
analysis of the proposed LS optimization is conducted, in 
which sufficient conditions and proofs for the convergence 
of the proposed algorithms are derived. The performance of 
the proposed scheme is assessed via simulations for MIMO 
equalization applications. The main contributions of this work 
are summarized as follows: 

1) A reduced-rank MIMO equalization scheme and a design 
approach for both decision feedback and linear structures; 

2) Reduced-rank LS expressions and recursive algorithms for 
parameter estimation; 

3) An algorithm for automatically adjusting the model order; 

4) Analysis and convergence proofs of the proposed algo- 
rithms. 

5) A study of MIMO reduced-rank equalization algorithms. 
This paper is structured as follows: The MIMO system 

and signal model is described in Section II. The proposed 



adaptive MIMO reduced-rank equalization structure is in- 
troduced along with the problem statement in Section III. 
Section IV is devoted to the development of the LS estimators, 
the computationally efficient RLS algorithms and the model 
order selection algorithms. Section V presents an analysis and 
proofs of convergence of the proposed algorithms. Section 
VI discusses the simulation results and Section VII gives the 
conclusions of this work. 

Notation: In this paper bold upper and lowercase letters 
represent matrices and vectors, respectively. (.)*, (.)*H, 
and (.) T shall represent complex conjugate, complex conjugate 
transpose (Hermitian), inverse and transpose, respectively. 
tr(.) is the trace operator of a matrix. Reduced-rank vectors 
and matrices are given with the addition of a bar (7) and 
estimated symbols are denoted by the addition of a hat (".). 

II. MIMO System and Signal Model 

In this section we present MIMO communications system 
and signal model and describe its main components. The 
model in this section is intended for describing a general 
MIMO system in multipath channels. However, it can also 
serve as a model for broadband MIMO communications 
systems with guard intervals including those based on or- 
thogonal frequency-division multiplexing (OFDM) ll37l . 11381 
and single-carried (SC) modulation with frequency-domain 
equalization 11391 . 

Consider a MIMO system with Nt antennas at the transmit- 
ter and Nr antennas at the receiver in a spatial multiplexing 
configuration, as shown in Fig. Q] The system is mathemati- 
cally equivalent to that in Jg). The signals are modulated and 
transmitted from Nt antennas over multipath channels whose 
propagation effects are modelled by finite impulse response 
(FIR) filters with L p coefficients, and are received by Nr 
antennas. We assume that the channel can vary during each 
packet transmission and the receiver is perfectly synchronized 
with the main propagation path. At the receiver, a MIMO 
equalizer is used to mitigate IAI and ISI and retrieve the 
transmitted signals. 
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Fig. 1. MIMO system model. 

The signals transmitted by the system at time instant i can 
be described by x[i] = \x\\i\ . . . X7v T [i]] T , where Xj[i], j = 
1 , ... Nt are independent and identically distributed symbols 
of unit variance. The demodulated signal received at the fcth 
antenna and time instant i after applying a filter matched to 
the signal waveform and sampling at symbol rate is expressed 
by 

N T Lp — 1 

Vk[i]='^2 X! ^j,kA i ] x A i ~^]+ n k[i}, f°r k=l,...,N R , 
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where hj t k,i[i] is the sampled impulse response between trans- 
mit antenna j and receive antenna k for path I, and n k [i] are 
samples of white Gaussian complex noise with zero mean and 
variance a 2 . By collecting the samples of the received signal 
and organizing them in a window of L symbols (L > L p ) for 
each antenna element, we obtain the LNr x 1 received vector 



y[i] = H[i]xT[i] + n[i 

T 



(2) 

where y[i] = [y'{[i] ... yjj contains the signals 

collected by the Nr antennas, the L x 1 vector y k [i] = 
[y k [i] ■ ■ ■ Vk[i - L + 1]] , for k = 1, . . . , N R , contains 
the signals collected by the fcth antenna and are organized 
into a vector. The window size L must be chosen according 
to the prior knowledge about the delay spread of the multipath 
channel J45]|. The LN R x LN T MIMO channel matrix H[i] 



H 2 ,i[i] 



H 2>2 [i] 



H2,N T [i] 



H 



H 



H 



N R ,N T \ l \ 



(3) 



where the L x L matrix Hj t k[i] are Toeplitz matrices with 
the channel gains organized in a channel vector hj^[i] = 
[hj,k,i[i] ■■■ hj.k.L -i] T that is shifted down by one position 
from left to right for each column, and which describes the 
multi-path channel from antenna j to antenna k. The elements 
hj,k,i[i], for 1 = 0, . . . , L p , of hj.k[i] are modelled as random 
variables and follow a specific propagation channel model 
B31 . as will be detailed in the Section VI. The LN T x 1 
vector Xx[i] = \p^\\f\ ■ ■ ■ ^a^H] i s composed by the data 
symbols transmitted from the Nt antennas at the transmitter 



with Xj [i] = [xj [i] 



i—L+1]] being the zth transmitted 



block with dimensions L x 1. The LNr x 1 vector n[i] 
is a complex Gaussian noise vector with zero mean and 
i?[n[i]n H [i]] = a 2 I, where (-) T and (-) H denote transpose 
and Hermitian transpose, respectively, and E[-] stands for 
expected value. 

III. Proposed Adaptive Reduced-Rank MIMO DFE 
and Problem Formulation 

We present the proposed reduced-rank MIMO equalization 
structure and state the main design problem of reduced- 
rank MIMO equalization structures. Both decision feedback 
(DF) and linear equalization structures can be devised by 
adjusting the dimensions of the filters and the use of feedback. 
We shall start with the description of the DF structure and 
then obtain the linear scheme as a particular case. In the 
proposed MIMO reduced-rank DF equalizer (DFE), the signal 
processing tasks are carried out in two stages, as illustrated 
in Fig. The proposed scheme employs two sets of filters 
and stacks the decision and the input data vectors for joint 
processing. The decision feedback strategy adopted in this 
work is the parallel scheme reported in (0, lTT3l . which firstly 
obtains the decision vector XTj[i) with linear equalization 
and then employs £C7Vj[i] to cancel the interference caused 
by the interfering streams. A decision delay <5d cc is assumed 
between the symbols transmitted and the x^j [i] obtained 
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after the decision block. The parallel strategy outperforms the 
successive one that uses a sequential procedure of equalization 
and interference cancellation Q, J8). 

Let us consider the design of the proposed MIMO reduced- 
rank equalizer using the structure shown in Fig. |2] The M x 1 
input data vector r[i] to the proposed equalizer is obtained by 
stacking the LNr x 1 received vector y[i] and the B(Nt — 
1) x 1 vector of decisions xj-jji] for stream j and is described 
by 

X T j[l\ 

LNr + B(Nt — 1) represents the number of 
samples for processing. The B(Nt — 1) x 1 vector of decisions 
xtj [i] = [xj [i] ... Xj [i — B + 1]] T for the jth stream takes 
into account B decision instants for the feedback and ex- 
cludes the jth detected symbol to avoid cancelling the desired 
symbol. The Nt x 1 vector of decisions is given by x[i] = 
[x\ [i] ... xn t [i]] T , whereas the Nt — lxl vector of decisions 
that excludes stream j and is employed to build xtj [i] is given 
by xj[i] = [xi[i] ... Xj-i[i] ■■■ x Nt [i\] t . 

Let us now consider an M x D transformation matrix 
SD,j[i] which carries out a dimensionality reduction on the 
received data rj[i] and shall exploit the low-rank nature of 
the data transmitted over stream j as follows 



where M 



r j \i] = SD, j [i\r j \i} > j = l,...,N 7 



(5) 



where D is the rank of the resulting equalization system. 

The resulting projected received vector fj[i] is the input 
to an estimator represented by the D x 1 vector Wj[i] = 
[Wji[i] Wj 2 [i] ... w{j£)[£]] T . According to the schematic 
shown in Fig. 12 the output of the proposed MIMO reduced- 
rank DFE is obtained by linearly combining the coefficients 
of Sr>,j[i] and Wj[i] for extracting the symbol transmitted 
from antenna j. Notice that all £>-dimensional quantities have 
a "bar". The proposed MIMO reduced-rank DFE output is 



>f[i]Si 



(6) 



From the outputs Zj[i] for j = 1, 2, . . . , Nt, we construct the 
vector z[i] = [zi[i] . . . Zj[i] . . . Zn t [i\\ . The initial decisions 
for each data stream are obtained without resorting to the 
feedback and are computed as follows 



vm 
o 



(7) 



where Q ( • ) represents a decision device suitable for the con- 
stellation of interest (BPSK, QPSK or QAM) and the vector of 
decisions is constructed as x[i] = [xi[i] ...&,-[£]... xn t 
and used to construct ir.jH and rj[i] as in (|4]i. The detected 



symbols x 



of the proposed reduced-rank MIMO DFE 



after the IAI and ISI cancellation are obtained by 

of 



x {f) \i] 



Q(z\i]) 



(8) 

The feedback employs B(Nt — 1) connections for cancelling 
the IAI and the other Nt — 1 data streams and the ISI from the 
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Fig. 2. Proposed MIMO reduced-rank decision feedback equalization structure. 



adjacent symbols. A reduced-rank MIMO linear equalizer is 
obtained by neglecting the feedback with decision processing 
of the structure in Fig. [2] 

The previous development suggests that the key aspect 
and problem to be solved in the design of reduced-rank 
MIMO equalization schemes is the cost-effective computation 
of the estimators Sd ,j [i] and ibj [i]. The transformation matrix 
SD,j[i] plays the most important role since it carries out the 
dimensionality reduction, which profoundly affects the perfor- 
mance of the remaining estimators and the MIMO equalizers. 
Methods based on the EVD flU, the MSWF (21 and the AVF 
l26l - ll28l were reported for the design of SD,j[i], however, 
they did not consider jointly the design of Sujii] and Wj[i] 
via alternating optimization recursions. In the next section, we 
present the reduced-rank least squares (LS) algorithms and 
their recursive versions for the design of Sr>,j[i] and Wj[i] 
used in the proposed MIMO equalization structure. 

IV. Proposed Reduced-Rank Least Squares Design 
and Adaptive Algorithms 

In this section, we present a joint iterative exponentially 
weighted reduced-rank LS estimator design of the parame- 
ters Sn,j[i] and Wj[i] of the proposed MIMO reduced-rank 
DFE. We then derive computationally efficient algorithms for 
computing the proposed LS estimator in a recursive way 
and automatically adjusting the model order. The deficient 
exchange of information between the dimensionality reduction 
task and the reduced-rank estimation verified in previously 
reported algorithms E1 - II281 is addressed by the alternated 
procedure that updates SD,j[i] and u>j-[i]. Specifically, the 
expression of Sd ,j [i] is a function of Wj [i] and vice versa, and 
this allows the coefficients to be computed via an alternating 
procedure with exchange of information in both ways (from 
SD,j[i] to Wj[i] and the other way around). Our studies and 
numerical results indicate that this approach is more effective 
than the MSWF |23) and the AVF EH) algorithms. In addition, 
the rank reduction is based on the joint and iterative LS 
minimization which has been found superior to the Krylov 
subspace, as evidenced in the numerical results. This allows 
the proposed method to outperform the MSWF and the AVF. 
We have opted for the use of one cycle (or iteration) per time 
instant in order to keep the complexity low. We also detail 
the computational complexity of the proposed and existing 
algorithms in terms of arithmetic operations. 



A. Reduced-Rank Least Squares Estimator Design 

In order to design Su.jii] and «Jj[£], we describe a joint 
iterative reduced-rank LS optimization algorithm. Consider 
the exponentially-weighted LS expressions for the estimators 
SD,j[i] and Wj[i) via the cost function 

i 

1=1 

(9) 

where < A < 1 is the forgetting factor. 

The proposed exponentially-weighted LS design corre- 
sponds to solving the following optimization problem 

In order to solve the problem in (fTOl . the proposed strategy is 
to fix a set of parameters, find the other set of parameters that 
minimize (O and alternate this procedure between the two sets 
SD,j[i] and ibj[i]. By minimizing © with respect to So,j[i], 
we obtain 

S Dtj [{\ = RJ 1 [i]PD,j H-Rfe,- [i - 1], (ID 

where the M x D matrix Pd j [i] = 

EL X^iljrMwfli ~ 1], = £l=i X~ l rj[l\rf[l], 
(■)' denotes the Moore-Penrose pseudo-inverse and the 
D x D matrix Rw 3 [i — 1] = Wj[i — l]wf[i — 1]. Since 
Rw [i — 1] is a rank-1 matrix, we need to either compute 
the pseudo-inverse or introduce a regularization term in the 
recursion FLw j [i — 1] = £; = J X l ~ l Wj[l]w^[l]. We have 
opted for using the latter with the initial regularization factor 
Rjjj. [0] = SI for numerical and simplicity reasons. 

By minimizing (O with respect to the reduced-rank 

estimator becomes 

w j \i}=Rj 1 \i]p j \i], (12) 

where p. [i] S%j [i] ZUi Pl»*i W 

Ya^^^WM], and the D x D reduced- 
rank correlation matrix is described by Rj[i] = 
S D,j[i\ E;=i ^~ l rj[l\rf[l\S Dlj [i]. 

The equation with the associated sum of error squares (SES) 
is obtained by substituting the expressions in (fTTT i and (fT2l into 
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the cost function (0, and is given by 

SES = al. -w?\i]S% Mp[i 



H r„-i e 



j L°J — ,j rjfrj P I*, 

w H [i]S" MRj[i\S D j[i\Wj[i 



(13) 



where cr^. = Ya=i \ x j ['] I 2 - Note that the expressions 
in ( fTTT > and ([T2l are not closed-form solutions for u>j[i] and 
SD.j[i] since they depend on each other and, thus, they have 
to be alternated with an initial guess to obtain a solution. The 
key strategy lies in the joint optimization of the estimators. 
The rank D must be set by the designer to ensure appropriate 
performance. The computational complexity of calculating 
( fTTT i and ( fT~2l > is cubic with the number of elements in the 
estimators, namely, M and D, respectively. In what follows, 
we introduce efficient RLS algorithms for computing the 
estimators with a quadratic cost. 

B. Reduced-Rank RLS Algorithms 

In this part, we present a recursive approach for efficiently 
computing the LS expressions developed in the previous sub- 
section. Specifically, we develop reduced-rank RLS algorithms 
for computing Wj[i] and Sflj[j]. Unlike conventional (full- 
rank) RLS algorithms that require the calculation of one 
estimator for the MIMO DFE, the proposed reduced-rank RLS 
technique jointly and iteratively computes the transformation 
matrix Sp and the reduced-rank estimator ti>j[i]. In order 
to start the derivation of the proposed algorithms, let us define 

i] = ^![<-i]> 



Pj[t\ 



(14) 



P Dij [i\ ± XP D ,j[i - 1] +x*\i]r j [i\w"[i - 1]. 
Rewriting the expression in (IPTl . we arrive at 

S DJ \i] = Rj 1 \i]P D , j \i]R£.[i-l] 

= P j [i\PD,MQ* ] li- 1 } 
= S D j[i - 1] + k 3 [i](x*[i]tf[i - 1] 
-r»\i}S Dij [i-l]), 

where the D x 1 vector tj [i — 1] = . [i — l]wj [i - 
M x 1 Kalman gain vector is 



kj \i] 



\- l Pj[i -l]r s 



l + X-WfWP^-l}^ 



(15) 



1], the 



(16) 



the update for the M x M matrix Pj [i] employs the matrix 
inversion lemma |fl9l 

Pj[t\ = X^Pjii - 1] - A-^j-Wrf [i]P 3 -[i - 1], (17) 

and the D x 1 vector tj[i — 1] is updated as 

A^QoJ* -!]«,•[*-!] 



tj[i-l) 



(18) 



The matrix inversion lemma is then used to update the D x D 
matrix Q^^ [i — 1] as described by 

g*,[i-l] = X-'Q^li^-X-H^-^wf [i-2]Q a [i-2]. 

(19) 



Equations (IT4T) -(fT9b constitute the part of the proposed 
reduced-rank RLS algorithms for computing Sojli]- 

In order to develop the second part of the algorithm 
that estimates tUj-[i], let us consider the expression in 
( fT2b with its associated quantities, i.e., the D x D matrix 
= Ei=i A i-i rj[Z]ff [I] and the D x 1 vector p^i] = 

Let us now define = Rj and rewrite pM\ as 

1] + x*[i]fj[i]. We can then rewrite (fl2l i as 



XPj [i 



follows 



w 3 \i] = *j[*W»] 



w, [i - 1] - k 3 \i\rf \i\wj + kj {i]x* [i] (20) 
w 3 [i - 1] + fcj [i] (xj [i] - rf [*]% [i - 1]) . 



By defining = Xj[i] — Wj[i — l]vj[i] we arrive at the 
proposed RLS algorithm for computing Wj [i] 



w A i ] = w j[i - 1] + fejMCj* W, 

where the D x 1 Kalman gain vector is given by 



l + A-iff [i]* J -[i-l]r i [i] : 



(21) 



(22) 



and the update for the matrix inverse <&[£] employs the matrix 
inversion lemma |[T9l 



= A-%[z - 1] - A-^.Wrf [z]*,^ - 1] 



(23) 



Equations d2TT >-(f23T> constitute the second part of the proposed 
algorithm that computes ibj [i] . The computational complexity 
of the proposed RLS algorithms is 0(D 2 ) for the estimation 
of Wj[i] and 0(M 2 ) for the estimation of Sdj^]- Since 
D « M for moderate to large L, Nr, Nt and B, as will 
be explained in the next section, the overall complexity is in 
the same order of the conventional full-rank RLS algorithm 

(0(m 2 )) Q3. 

C. Model-Order Selection Algorithm 

The performance of the LS and RLS algorithms described 
in the previous subsection depends on the model order or 
the rank D. This motivates the development of methods to 
automatically adjust D using an LS cost function as a mech- 
anism to control the selection. Prior methods for model order 
selection which use MSWF-based algorithms [23] or AVF- 
based recursions 11281 have considered projection techniques 
ll23l and cross-validation l28l approaches. Here, we focus 
on an approach that jointly determines D based on an LS 
criterion computed by the estimators Srjj[i] and Wj[i], where 
the superscript D denotes the rank used for the adaptation. The 
methods considered here (the proposed and existing ones 11231 . 
1 28 1) are the most suitable for model-order adaptation in time- 
varying channels. Other techniques such as the Akaike infor- 
mation criterion-based and the minimum description length 
do not lend themselves to time-varying situations and are 
computationally complex lfT9l . 

The key quantities to be updated are the transformation ma- 
trix Srijli], the reduced-rank estimator tUj-[i], and the inverse 
of the reduced-rank covariance matrix Pj [i] (for the proposed 
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RLS algorithm). Specifically, we allow the dimensions of 
SD,j[i) an d to vary from !),„;„ and D max , which are 

the minimum and maximum ranks allowed, respectively. It is 
important to note that only one recursion to obtain Pj[i] is 
computed with Z? max in order to keep the complexity low. 
Once Pj [i] is obtained, we perform a search for the best D 
for SD.j[i] and Wj[i] that require sub-matrices of Pj[i] for 
their computation. The transformation matrix So.jii] and the 
reduced-rank estimator w j [i] employed with this algorithm are 
illustrated by 



TABLE I 

Computational complexity of algorithms. 



and 



SMS,j[i] 



w D,j[i] = [ wij[i] w 2 ,j[i] ... w Dmin j[i] 

(24) 

The method for automatically selecting D of the algorithm is 
based on the exponentially weighted a posteriori least-squares 
type cost function: 

i 

CjiSoj [i\,w dJ [*])=£ W - <j W S d,i Mr, HI 2 - 

1=1 

(25) 

For each time interval i, we select the rank -D pt[«] which 
minimizes Cj(SD,j[i],WD,j[i]) an d the exponential weighting 
factor A is required as the optimal rank varies as a function 
of the data record. The transformation matrix So and the 
reduced-rank estimator w<i,j[i] are updated along with P[i] 
for the maximum allowed rank £> max and then the proposed 
rank adaptation algorithm determines the the best model order 
for each time instant i using the cost function in (f25|. The 
proposed model-order selection algorithm is given by 



D 



Jioptl 



arg 



mm 

Qin <d<r>„ 



C^Sd^ilwdjli}), (26) 



where d is an integer, D B 
maximum ranks allowed for the estimators, respectively. A 
small rank may provide faster adaptation during the initial 
stages of the estimation procedure, whereas a large rank 
usually yields a better steady-state performance. Our studies 
indicate that the range for which the rank D of the proposed 
algorithms have a positive impact on the performance of the 
algorithms is limited. Specifically, we have found that even for 
large systems (N R = Nt = 20,30,40,50,60) the rank does 
not scale with the system size and remains small. The typical 
range of values remains between D mm = 3 and D max = 8 for 
the system sizes examined (N R = N T = 20,30,40,50,60). 
This is an important aspect of the proposed algorithms because 
it keeps the complexity low (comparable to a standard RLS 
algorithm). For the scenarios considered in what follows, we 
set -D m in = 3 and Z) max = 8. In the simulations section, 
we will illustrate how the proposed model-order selection 
algorithm performs. 

D. Computational Complexity 

In this subsection, we illustrate the computational com- 
plexity requirements of the proposed RLS algorithms and 



A 1 ci/~\t"i t ri m 
rtl^UIlLIHIl 


rtUUlLlUlla 


IVlUlLipilL-clLIUIla 


Full-rank N91 

■ uii i aim | i ' | 


2M 2 


3M 2 




+M + 1 


+5M 


Pronoscd 


2M' 2 


3M 2 




-A/ + 4D 2 


+3M + 6D 2 




+MD +D+3 


+MD + 8D 


MSWF |23] 


DM 2, 


DM' 2 




+6D 2 -8D + 2 


+2DM + 3D 




+M 2 


+M 2 + 2 


AVF (28) 


DM 2 + 2M - 1 


4DM 2 




+5D(M - 1) + 1 


+4DM + AM 




+3(DM - l) 2 


+4D + 2 



compare them with those of existing algorithms. We also 
provide the computation complexity of the proposed and 
existing model-order selection algorithms. The computational 
complexity of the algorithms is expressed in terms of additions 
and multiplications, as depicted in Table I. For the proposed 
reduced-rank RLS algorithm the complexity is quadratic with 
M = LNr + B(Nt — 1) and D. This amounts to a complexity 
slightly higher than that observed for the full-rank RLS 
algorithm, provided D is significantly smaller than M, and 
significantly less than the cost of the MSWF-RLS (23l and 
the AVF |[28 1 algorithms. The complexity of the proposed 
model-order selection algorithm is given in Table II. 

In order to illustrate the main trends and requirements in 
terms of complexity of the proposed and existing algorithms, 
we show in Fig. [3] the complexity against the number of input 
samples M for the parameters D = 5, Nt = Nr, L = 8 
and B = 2. The curves indicate that the proposed reduced- 
rank RLS algorithm has a complexity significantly lower than 
the MSWF-RLS algorithm J23) and the AVF (28), whereas it 
remains at the same level of the full-rank RLS algorithm. 




Fig. 3. Complexity in terms of multiplications against number of 
input samples (M) with D = 5, N T = N R , L = 8, B = 2 . 



The computational complexity of the model-order selection 
algorithms including the proposed and the existing techniques 
is shown in Table II. We can notice that the proposed model- 
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TABLE II 

Computational complexity of model-order selection 
algorithms. 



The proposed LS and RLS algorithms can be stated as an 
alternating minimization strategy based on the sum of error 



Algorithm 

Proposed 


Additions 

2(Anax - Anin) + 1 


Multiplications 


Projection with 
Stopping Rule 1231 


2(2M - l)x 

(Anax — -^min) H~~ 1 


(M a + M + l)x 

(^max Dmin ~\~ 1) 


cv no 


(2M - l)x 

(2(Anax ~ Anin) + 1) 


(Anax - Anin + 1) X 
M+l 



G arg 



min SES(S'^ pt 7 ,w J [i]) 



%[i]earg min SES(S D , 3 [i},w° pt ), 

w° pt £w t [i] 



(28) 
(29) 



order selection algorithm is significantly less complex than the 
existing methods based on projection with stopping rule 11231 
and the CV approach l28l . Specifically, the proposed algorithm 
that uses extended filters only requires 2(£> max — -Dmin) 
additions, as depicted in the first row of Table II. To this 
cost we must add the operations required by the proposed 
RLS algorithm, whose complexity is shown in the second 
row of Table I using D max . The complexities of the MSWF 
and the AVF algorithms are detailed in the third and fourth 
rows of Table I. For their operation with model-order selection 
algorithms, a designer must add their complexities in Table I 
to the complexity of the model-order selection algorithms of 
interest in Table II. 

V. Analysis of The Proposed Algorithms 

In this section, we conduct an analysis of the proposed 
algorithms that compute the estimators SD,j[i) and Wj[i) 
of the proposed scheme. We first highlight the alternating 
optimization nature of the proposed algorithms and make use 
of recent convergence results for this class of algorithms fiOl . 
flJTJ. In particular, we present a set of sufficient conditions 
under which the proposed algorithms converge to the optimal 
estimators. This is corroborated by our numerical studies that 
verify that the method is insensitive to different initializations 
(except for the case when So j[i] is a null matrix which 
annihilates the received signal) and that it converges to the 
same point of minimum. We establish the global convergence 
of the proposed algorithm via induction and show that that the 
sequence of estimators SDj[i] an d w j[i] produces a sequence 
of outputs that is bounded and converges to the reduced-rank 
Wiener filter SEE]. 

A. Sufficient Conditions for Convergence 

In order to develop the analysis and proofs, we need to 
define a metric space and the Hausdorff distance that will be 
used extensively. A metric space is an ordered pair (Ai,d) 
where Ai is a non-empty set and d is a metric on At, i.e., a 
function d:MxM->R such that for any x, y, z, and Ai 
we have: 

a) d(x,y) > 0. 

b) d(x,y) = iff x = y. 

c) d(x,y) = d(y,x). 

d) d(x, z) < d(x, y) + d(y, z (triangle inequality). 

The Hausdorff distance measures how far two subsets of a 
metric space are from each other and is defined by 

dn{X,Y) = maxjsup inf d(x,y),swp inf d(x,y)} (27) 



where • and Wj correspond to the optimal values of 
Spjli] and respectively, and the sequences of compact 

sets {S D j[i]}i> and {w^i}}^ converge to the sets S_ D opt 
and Wj t , respectively. 

Although we are not given the sets S_ D t and w_- t di- 
rectly, we observe the sequence of compact sets {S_ D j[i]}i>o 
and {Wj [i]}i>o- The goal of the proposed algorithms is to find 
a sequence of Sujii] and Wj[i] such that 

lim SES(Su *[»],«><[«]) = SES(S'° pt ,,io° pt ) (30) 

In order to present a set of sufficient conditions under which 
the proposed algorithms converge, we need to so-called 
"three-point" and "four-point" properties iffOl . iBTI . Let us 
assume that there is a function / : A4 x M — > M such that 
the following conditions are satisfied: 

1) Three-point property (S"^, So,j, w° pt ): for all i > 1, 



■ opt 

D,j c ii_D.jL'-J' 

S Dtj G argmin ii) opt g 



™ 7 - [*] 



SES(^ p *,^ opt ) 



f{S°g;,S Dd ) + SES(5 D , J ,«;° pt ) < SESCS^.,^). 



2) Four-point property (S opt 



all i > 1, 



oopt 



w 



opt 



>D,j, 



Sd,j G SjjM, 



(31) 
w° pt ): for 



w° pt G w_j[i], and 
w D,j G argmin^opt^^ SES(5 D , 3 ■,, w° pt ) 

SES(Sg?., Wj ) < SES(S op ^w° pt ) + f(S op ^ S D ,j). (32) 

Theorem: Let {(S_D,j[i\iSLj[i])}i>o> Sjdj> liLj Pt be com - 
pact subsets of the compact metric space (Ai, d) such that 



qOpt 



w 



o p t 



(33) 



and let SES : M x M — > K be a continuous function. 

Now let conditions 1) and 2) hold. Then, for the proposed 
algorithms we have 



(34) 



lim SES(5 Bj [i],^[z]) = SES(S3 > ' tcf *) 
A general proof of this theorem is detailed in ROl . PT1 . 



B. Convergence to the Optimal Reduced-Rank Estimator 

In this subsection, we show that the proposed reduced- 
rank algorithm converges globally and exponentially to the 
optimal reduced-rank estimator ||20l , ll2TI . We assume that 
1 > A ^> (equal or close to one) , the desired product of the 
optimal solutions, i.e., w° pt = S C p t jW° p is known and given 

by Rf /2 [i](Rf /2 \i]pj[i]) VD HEED, where #7 1/2 H is 

the square root of the input data covariance matrix and the 
subscript 1 : D denotes truncation of the subspace. 
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In order to proceed with our proof, let us rewrite the 
expressions in ( fTTT > and ( TT2T > for time instant as follows 

ft, (KS/,.., U 7?„. [0] = P DJ [0] = Pj [0]wf[0], (35) 

R 3 [OjiVj [1] = sZj [0]Rj [0)S D ,j [0]wj [1] - v 3 [0] , (36) 

Using ( f35l > we can obtain the following relation 



sZMRMpMwfM, 



At this point, we need to establish that the norm of So.j [i] for 
all i is both lower and upper bounded, i.e., < ||S , £>. :) [i]|| < 
oo for all i, and that Wj[i] = Sjj t j[i]wj[i] approaches 
w j,opt[i] exponentially as i increases. Due to the linear map- 
ping, the boundedness of Sdj [i] is equivalent to that of Wj [i]. 
Therefore, we have upon convergence w^[i]wj[i — 1] = 
wj[i — — 1]. Since Hio^fiJiUj-^ — 1]|| < \\wj[i — 

l]\\\\wj[i]\\ and — 



(37) ^ e re l at i° n w f\f\ w i\?' ~ 1] = w f[i ~ ^\ w j[i — 1] implies 



Substituting the above result for R Wj [0] into the expression in 
fl35l l we get a recursive expression for Sd,j[0] 

(38) 



and hence 

[oo] 1 1 > > \\wj[ 



(47) 



{S^R^OjSnM) \ 



Using fl36l ) we can express Wj [1] as 



SnJObJO], (39) 



In order to show that the upper bound ||u> 3 -[oo]|| is finite, 
let us express the M x M matrix Qj[i] as a function of the 

w jt2 [i] 



M x 1 vector Wj[i] 
Aj,i 



and the M x M matrix 



A = 
of w 



Employing the relation Wj[l] = S D.j[]]wj[l], we obtain 

U7 3 -[1] = R.iirM^fmsl^RA^P^wfll})- 1 - 

■{Sl^R^So^y^Sl^RMSDAO^sZm^ 

(40) 

More generally, we can express the proposed reduced-rank LS 
algorithm by the following recursion 



Substituting the previous expressions 



t] and Aj into Q ■ [i] given in (1461) . we obtain 



A" «,-![()] 

trfalOlA?'- 2 ^^)- 1 



[i]-^ [*]*&f [i] [»] R 3 [i\ Pj [i]wf W) 

- 1 e ff 



Using the matrix identity (A + B) 







<i[°]A?r a 


, (48) 


< 2 [o]a?; 2 - 2 




= A" 1 - A" 





l A 1 to the decomposed QJi] in d48| ) and making 



(Sg^i - l]J^[t - l]5z,j[i - I]) S% tj [i - l] Pj [i - 1] 

(41) 



i large, we get 

Q 3 \i] =diag (!_._! (L.jT 



O(eW). 



(49) 



M-D 



Since the optimal reduced-rank filter can be described by the 
EVD of JJT 1/2 [i] Pj .[i] |2Q), ED, where B7 1/2 [i] is the square 
root of the covariance matrix and Pj-[i] is the cross- 

correlation vector, then we have 



where e[i] = (A r +i/A r ) with A r +i and A r are the (r + 
l)th, the rth largest singular values of R- 1 [i]pM] and O(-) 



H7 1/2 [i]p,-H 



^A^f^-W, 



(42) 



denotes the order of the argument. From (1491 . it follows that 
for some positive constant k, we have ||iOj[i]|| < 1 
From d43T >. we obtain 



ke\i] 



where Aj is a M x M diagonal matrix with the eigenvalues of 
Rj and <I>j is a M x M unitary matrix with the eigenvectors 
of Rj. Assuming that there exists some Wj[0] such that the 
randomly selected SD,j[0] can be written as fl2~Tl 



s Dtj [o]=R; 1/2 mj™M- 



(43) 



Using d43l and ( l42b in ( |4TI ). and manipulating the algebraic 
expressions, we can express (|4"TT i in a more compact way that 
is suitable for analysis, as given by 

Wj[i] = A*u>j[i-l](u>f [i-ilAjWjli-l])' 1 ™? [i-l]wj[i-l\. 

(44) 

The above expression can be decomposed as follows 



< 



Q i [oo]||...||Q J .[2]|| \\Q 3 \ 

t»i[0]|in(l + Ae[i]) 
i=l 



\Qi 



w 



exp {^2 l °9( i + fce W) 

i=l 



(50) 



| cxp 



| exp 



(E* e w) 



1 - (A r+ l/A r ) : 



i» J [i] = Q i [i] Qj[i-1] ■■■ QM 



where 



Q,.[i] = Af ^.[0]«[0]Af- 2 ^[0])- 1 <[0]A 



, 2i— 2 



(45) With the development above, the norm of t«j [i] is proven to 
be both lower and upper bounded. Once this is established, 
the expression in (l4T~b converges for a sufficiently large i to 

(46) the reduced-rank Wiener filter. This is verified by equating the 



9 



terms of (l44l which yields 

Wj \i] = Rj [i] " V,- [i]wf [i] (Sg^ [i]p 3 - [i]t»f [i]) _1 ■ 



N T =4, N R =8, L=8, SNR=15 dB, f T=0.0001, Data record = 500 symbols 



(S^[i - ljRAi - l]S Dd [i - l\) S^[i - l] Pj [i - 



R 



-1/2 



i]*i,i A i,i*fiPi[*1 + 0(e[i]), 



(51) 



where <&i is a A/ x D matrix with the D largest eigenvectors of 
Rj[i] and A 3i i is a D x D matrix with the largest eigenvalues 
of Rj[{\. 

VI. Simulation Results 

In this section, we evaluate the bit error rate (BER) 
performance of the proposed MIMO equalization structure, 
algorithms and existing techniques, namely, the full-rank (9), 
the reduced-rank MSWF J23], and the AVF J28) techniques 
for the design of the receivers. For all simulations and the 
proposed reduced-rank RLS algorithm, we use the initial val- 
ues u^O] = [1,0,..., 0] and S DJ [0] = [I D Dx(M -d)) T - 
For the next experiments, we adopt an observation window 
of L = 8, the multipath channels ( the channel vectors 
hj,k[i] = [hj,k,i[i] ••• %fc,L p -i] T ) are modelled by FIR 
filters with L p coefficients spaced by one symbol and the 
system employs QPSK modulation. The channel is time- 
varying over the transmitted packets, the profile follows the 
UMTS Vehicular A channel model 04) with L p = 5 and 



the fading is given by Clarke's model H31 . We average the 
experiments over 1000 runs and define the signal-to-noise ratio 
(SNR) as SNR = 101og 10 ^ji, where a 2 x is the variance of 
the transmitted symbols and a 2 is the noise variance. The 
adaptive MIMO equalizers employ N T = 4, B = 4, L = 8 
and Nji = 8 in a spatial multiplexing configuration, leading to 
estimators at the receiver with M = LNr + B{Nt — 1) = 72 
coefficients. The adaptive RLS estimators of all methods are 
trained with 250 symbols, employ A = 0.998 unless otherwise 
specified, and are then switched to decision-directed mode. 

A. Convergence Performance and Impact of Rank 

In the first experiment, we consider the BER performance 
versus the rank D with optimized parameters (forgetting factor 
A = 0.998) for linear MIMO equalizers. The curves in Fig. 4 
show that the best rank for the proposed scheme is D = 4, 
which is the closest among the analyzed algorithms to the 
optimal linear MMSE that assumes the knowledge of the 
channel and the noise variance. In addition, it should be 
remarked that our studies with systems with different sizes 
show that the optimal rank D does not vary significantly 
with the system size. It remains in a small range of values, 
which brings considerably faster convergence speed. However, 
It should also be remarked that the optimal rank D depends 
on the data record size and other parameters of the systems. 

The BER convergence performance versus the number of 
received symbols for MIMO decision feedback equalizers with 
optimized but fixed ranks is shown in Fig. 5. The results 
show that the proposed scheme has a significantly faster 
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Fig. 4. BER performance versus rank (D) for linear MIMO equalizers. 
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Fig. 5. BER performance versus number of received symbols. 



convergence performance and obtains good gains over the best 
known schemes. The plots show that the proposed reduced- 
rank MIMO equalizer extends the dimensionality reduction 
and its benefits such as fast convergence and robustness against 
errors to the MIMO equalization task. The proposed RLS 
estimation algorithm has the best performance and is followed 
by the AVF, the MSWF, and the full-rank estimators. Note 
that the BER of the considered techniques will converge to the 
same values if the number of received symbols is very large 
and the channel is static. 

B. Performance with Model Order Selection 

As previously mentioned in Section IV, it is possible to 
further increase the convergence speed and enhance the 
tracking performance of the reduced-rank algorithms using 
an automatic model-order selection algorithm. In the next 
experiment, we consider the proposed reduced-rank structures 
and algorithms with linear and DF equalizers and compare 
their performance with fixed ranks and the proposed automatic 
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model-order selection algorithm developed in Section IV.C. 
The results illustrated in Fig. 6 show that the proposed 
model-order selection algorithm can effectively speed up the 
convergence of the proposed reduced-rank RLS algorithm and 
ensure that it obtains an excellent tracking performance. 
In what follows, we will consider the proposed model-order 
selection algorithm in conjunction with the proposed reduced- 
rank RLS algorithm, and for a fair comparison we will equip 
the MSWF and the AVF algorithms with the rank adaptation 
techniques reported in ||231 and J28), respectively. 



N T =4, N R =8, L=8, B=4, SNR = 1 5 dB, f d T=0.0001 

10° □ 1 i 




Number of received symbols 



Fig. 6. BER performance versus number of received symbols for proposed es- 
timation algorithms and structures. The performance of the proposed reduced- 
rank algorithms is shown for the proposed model-order selection algorithm 
(solid lines), for D = 3 (dotted lines) and for D = 8 (dash-dotted lines). 



C. Performance for Various SNR and f^T Values 

The BER performance versus the signal-to-noise ratio 
(SNR) for MIMO decision feedback equalizers operating with 
the automatic model-order selection algorithms is shown in 
Fig. 7. The curves show a significant advantage of reduced- 
rank algorithms over the full-rank RLS algorithm. Specifically, 
the reduced-rank AVF and MSWF techniques obtain gains 
of up to 3 dB in SNR for the same BER over the full-rank 
algorithm, whereas the proposed reduced-rank RLS algorithm 
achieves a gain of up to 3 dB over the AVF, the second best 
reduced-rank algorithm. The main reasons for the differences 
in diversity order are the speed and the level of accuracy of the 
parameter estimation of the proposed and existing methods. If 
we increase the number of received symbols to a very large 
value then the diversity order attained by the different analyzed 
algorithms would be the same as verified in our studies. 

In order to assess the performance of the reduced-rank al- 
gorithms for different fading rates, we consider an experiment 
where we measure the BER of the proposed and analyzed 
algorithms against the normalized fading rate f^T in cycles 
per symbol, where is the maximum Doppler frequency and 
T is the symbol rate. It should be noted that the forgetting 
factor A was optimized for each value of f^T in this 
experiment. In practice, a designer could employ a mechanism 



N T =4, N R =8, L=8, B=4, f d T=0.0001 , Data record=1500 symbols 
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Fig. 7. BER performance versus SNR . 



to automatically adjust A. The results of this experiment are 
shown in Fig. 8, where the advantages of the reduced-rank 
algorithms and their superior performance in time-varying 
scenarios is verified again. 



N T =4, N R =8, L=8, B=4, SNR = 15 dB, Data record = 1500 symbols 
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Fig. 8. BER performance versus the normalized fading rate f^T . 



D. Performance for MIMO-OFDM Systems 

In the previous experiments, we considered the proposed 
MIMO equalization structure and algorithms for time-varying 
channels that dynamically change within a packet transmis- 
sion, thereby requiring the adaptive equalization techniques so 
far described. At this point, it would be important to address 
two additional issues. The first is to account for the gains 
of reduced-rank techniques over full-rank methods when the 
order of the estimators changes. Another important aspect to 
be investigated is the applicability of the proposed reduced- 
rank techniques to broadband communications such as MIMO- 
OFDM systems EH, E3- Even though in MIMO-OFDM 
systems the frequency selective channels are transformed into 
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frequency flat channels, there is still the need to perform 
spatial equalization. We consider an experiment with a MIMO- 
OFDM system in which the data streams per sub-carrier are 
separated by MIMO linear equalizers equipped with full-rank 
and reduced-rank algorithms and the channels change at each 
OFDM block. The system has = 64 sub-carriers and 
employs a cyclic prefix that corresponds to C = 8 symbols. 
The channel profile is identical to the model employed for the 
previous experiments and the fading is independent for each 
stream. The Nr x 1 received data vector for the ?ith subcarrier 
is given by 

r n [i] = H n [i]x n [i] +n n [i], n = 1, 2, . . . , N, (52) 

where the Nr x Nt channel matrix H „ [i] contains the channel 
frequency response gains at the nth tone, the Nt x 1 data 
vector x n [i] corresponds to the symbols transmitted by the 
Nt antennas over the nth subcarrier and the Nr x 1 vector 
n n [i) represents the noise vector at the nth tone. 

We employ the proposed MIMO linear equalization scheme 
for spatial equalization on a per subcarrier-basis fl37l for the 
OFDM symbols with the proposed and analyzed reduced-rank 
estimation algorithms. The BER is plotted against the number 
of antennas in a MIMO-OFDM system that has N R = N T . 
The results in Fig. 9 show that the advantages of reduced- 
rank algorithms are more pronounced for larger systems, in 
which the training requirements are more demanding in term 
of training data for full-rank RLS algorithms. 



SNR = 15 dB, f T=0.0001, Data record = 1500 symbols, Training = 250 symbols 




Fig. 9. BER performance versus the number of antennas . 

The advantages of the reduced-rank estimators are due to the 
reduced amount of training and the relatively short data record 
(packet size). Therefore, for packets with relatively small 
size, the faster training of reduced-rank LS estimators will 
lead to superior BER to conventional full-rank LS estimators. 
As the length of the packets is increased, the advantages of 
reduced-rank estimators become less pronounced for training 
purposes and so become the BER advantages over full-rank 
estimators. In comparison with the MSWF and AVF reduced- 
rank schemes, the proposed scheme exploits the joint and 
iterative exchange of information between the transformation 



matrix and the reduced-rank estimators, which leads to bet- 
ter performance. The gains of the reduced-rank techniques 
over full-rank methods for MIMO-OFDM systems are less 
pronounced than those observed for narrowband MIMO sys- 
tems with multipath channels. This is because the number 
of coefficients for estimation is significantly reduced. If we 
increase the number of antennas in MIMO-OFDM systems 
to a large value, then the gains of reduced-rank techniques 
become larger. 

VII. Concluding Remarks 

This paper has presented a study of reduced-rank equal- 
ization algorithms for MIMO systems. We have proposed 
an adaptive reduced-rank MIMO equalization scheme and 
algorithms based on joint iterative optimization of adaptive 
estimators. We have developed LS expressions and efficient 
RLS algorithms for the design of the proposed reduced-rank 
MIMO equalizers. A model-order selection algorithm for auto- 
matically adjusting the model order of the proposed algorithm 
has also been developed. An analysis of the convergence of 
the proposed algorithm has been carried out and proofs of 
global convergence of the algorithms have been established. 
Simulations for MIMO equalization applications have shown 
that the proposed schemes outperforms the state-of-the-art 
reduced-rank and the conventional estimation algorithms at 
a comparable computational complexity. Future work and 
extensions of the proposed scheme may consider strategies 
with iterative processing |[L3l . (25], ll46) with the aid of 
convolutional, Turbo and LDPC codes, detection structures 
which attain a higher diversity order ifPJI . lTT~4l and their 
theoretical analysis. 
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